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Abstract 

This paper deals with redundancy allocation problem of 
complex system with ten subsystems connected in parallel 
with identical components. The reliability of each 
component is considered to be imprecise. This impreciseness 
has been represented by different fuzzy numbers. Then the 
corresponding problem has been transformed to crisp 
constrained optimization problem using different 
defuzzification methods. The transformed problem has been 
converted into unconstrained optimization problem by Big- 
M penalty technique, then solved by real coded genetic 
algorithm with advanced operators. Finally, to illustrate the 
methodology, a numerical example has been applied for 
different defuzzification methods and the obtained results 
have been compared. 
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Introduction 

Due to the development of modern technology, the 
design of a system is dependent on the selection 
of components and configurations to meet the 
requirement as well as the performance. For a system 
with known cost, reliability, weight, volume and other 
parameters, the corresponding design problem is 
known as combinatorial optimization problem. The 
well known reliability design problem is referred to as 
redundancy allocation problem whose fundamental 
objective is to determine the number of redundant 
components that maximize the system reliability or 
minimize the system cost subject to the several 
resource constraints. As redundancy allocation problem 
is a nonlinear integer programming problem, it can 



not be solved by direct/indirect or mixed search 
methods due to discrete search space. This redundancy 
allocation problem is NP-hard and has been well 
discussed by Tillman et al. (1977a) and Kuo and 
Prasad (2001). In the existing literature, it is observed 
that several methods, viz. heuristic methods 
(Nakagawa and Nakashima 1977; Kim and Yum 1993; 
Kuo et al. 1978; Aggarwal and Gupta 2005; Ha and 
Kuo 2006), branch and bound method (Kuo et al. 1987; 
Tillman et al. 1977; Sun and Li 2002 ; Sung and Cho 
1999), reduced gradient method (Hwang et al. 1979), 
integer programming method (Misra and Sharma 
1991), dynamic programming (Nakagawa et al. 1981; 
Hikita et al. 1992) etc. were used to solve such 
redundancy allocation problem. However, these methods 
have some advantages as well as disadvantages. 
Dynamic programming technique can not be used in 
solving complicated problem as it can not be 
decomposable according to 'Bellman's optimality 
principle'. In branch and bound method, the 
effectiveness depends on the sharpness of the bound 
and required memory increases exponentially with the 
problem size. As a result, with the development of 
genetic algorithm and other evolutionary algorithms, 
most of the researchers are motivated to use these 
methods in solving redundancy allocation problem. 
These methods provide more flexibility, require less 
assumptions on the objective as well as constraints 
and more efficient irrespective of whether the search 
space is discrete or not. 

In almost all the studies mentioned earlier, the 
parameters of redundancy allocation problem have 
usually been taken to be precise values. This means that 
complete probabilistic information about the system is 
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known. In this situation, it is commonly assumed that 
the value of probability is perfectly determinable for 
every event. However, in reality it is not true. There are 
not sufficient data available in most of the cases when 
the system is either new or it exists only as a project. In 
this case the data/ information can not be collected 
precisely due to human errors, improper storage 
facilities and other unexpected factors related to 
environment. So the reliability of a component of a 
system will be an imprecise number. To take the 
problem with such imprecise numbers, generally 
stochastic, fuzzy and fuzzy-stochastic approaches are 
applied and the corresponding problems are converted 
to deterministic problems for solving them. 

In this paper, redundancy allocation problem of 
complicated system with ten subsystems has taken 
into consideration. Each subsystem is connected 
parallelly with identical components. The reliability of 
each component is considered to be fuzzy number. 
Then the corresponding problem has been transformed 
to crisp problem using different defuzzification 
methods. The transformed problem has been converted 
into unconstrained optimization problem by Big-M 
penalty technique, then solved by real coded genetic 
algorithm for integer variables with tournament selection, 
intermediate crossover and one neighbourhood 
mutation. Finally, to illustrate the methodology, a 
numerical example has been listed for different 
defuzzification methods and the computed results 
have been compared. 

Representation of Fuzzy Number 

The word 'fuzzy' was first introduced by Zadeh in 
1965 in his famous research paper "Fuzzy Sets" as a 
mathematical way to represent impreciseness/ 
fuzziness or vagueness. The approach of fuzzy set is 
an extension of classical set theory and used in fuzzy 
logic. In classical set theory, the membership of each 
element in relation to a set is assessed in binary terms 
according to a crisp condition; and an element either 
belongs to or not to the set. By contrast, a fuzzy set 
theory permits the gradual assessment of the 
membership of each element in relation to a set; which 
is discussed with the aid of a membership function. 
Fuzzy set is an extension of classical set theory since, 
for a certain universe, a membership function may act 
as an indicator function, mapping all elements to 
either 1 or 0, as in the classical notation. He used this 
word to generalize the mathematical concept of the set 
to one of fuzzy set or fuzzy subset, where in a fuzzy 
set, a membership function is defined for each element 



of the referential set. 

Fuzzy Set: A fuzzy set A in a universe of discourse X 
is defined as the set of pairs 

A = {(x,jU A (x)):xeX], 
where n~ k : X -> [0,1] is a mapping and n~ k (x) is called 

the membership function of A or grade of 
membership of x in A . 

Convex Fuzzy Set: A fuzzy set A is called convex if 
and only if for all 

Xj,x 2 e X, jU^(Ax l +(l-A)x 2 ) > min{//^(x 1 ), / u^(x 2 )}, 

where /Le[0,l]. 

Support of a fuzzy set: The support of fuzzy set A 
denoted by S(A) is the crisp set of all x e X such that, 
ju A (x)>0. 

a -level Set: The set of elements that belong to the 
fuzzy set A at least to the degree a is called the a - 
level set or a -cut and is given by 

A a ={xeX:ju A (x)>a}. 

If A a = {x e X : jU A (x) > a} , it is called strong a -level 
set or strong a -cut. 

Normal Fuzzy Set: A fuzzy set A is called a normal 
fuzzy set if there has at least one x e X such 
that n~ A (x) = 1 . 

Fuzzy Number: A fuzzy number is a fuzzy set which 
is both convex and normal as well as a special case of a 
fuzzy set. Different definitions and properties of fuzzy 
numbers are encountered in the literature but they all 
agree on that a fuzzy number represents the notion of 
a set of real numbers 'closer to a ' where ' a ' is the 
number being fuzzified. 

Triangular Fuzzy Number (TFN) 

A triangular fuzzy number A is represented by the 
triplet (a 1 ,a 2 ,a 3 ) and defined by its continuous 
membership function /j^(x) : X -> [0,1] given by 
x-a, 

it Oj < x < a 2 

if x = a 2 
if a 2 ^ x < a 3 
otherwise 



fl 2 - flj 

1 



a-, -a 



Parabolic Fuzzy Number (PFN) 

A parabolic fuzzy number A is represented by the 
triplet (a 1 ,fl 2 ,fl 3 ) and defined by its continuous 
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membership function ju^(x) : X — > [0,1] given by 



Mi ( x ) : 



v"2 "iy 



x-a 



if flj < x < a 2 



if x = cu 



if a 2 < x < fl 3 



a 3 -a 2 ^ 
otherwise 
where, L(x) and /?(x) are continuous functions of x. 
Moreover, L(x) is strictly monotonic increasing and 
R(x) strictly monotonic decreasing function of x in 
flj < x < a 2 and a 3 < x < a 4 respectively. 

Trapezoidal Fuzzy Number (TrFN) 

A trapezoidal fuzzy number A is represented by the 
quadruplet (a[,a 2 ,a 3 ,a 4 ) an d defined by its continuous 
membership function ju^(x) : X -> [0,1] given by 
x — cu 

— it < x < a 2 



Mi W = 



if a 2 < x < a 3 



if a 3 < x <a 4 



otherwise 



Generalized Fuzzy Number 

The generalized fuzzy number A with membership 
function jU^(x) (Fig. 1) exhibits a fuzzy subset of the 

real line □ , where 



Ml (*) : 



L(x) flj < x < a 2 



1 



a 2 < x < a 3 



R(x) a 3 < x <a 4 
otherwise 




FIG. 1 GENERALIZED FUZZY NUMBER 

Defuzzification 

Defuzzification is the process of producing a 
quantifiable result in fuzzy logic, given the fuzzy sets 
and the corresponding degrees of membership. There 



are several defuzzification techniques available in the 
existing literature. However, the common and useful 
ones are as follows: 

(i) Centre of Area (COA) 

(ii) Bisector of Area (BOA) 

(iii) Smallest of Maxima (SOM) 

(iv) Largest of Maxima (LOM) 

(v) Mean of Maxima (MOM) 

(vi) Regular Weighted Point (RWP) 

(vii) Graded Mean Integration Value (GMIV) 

(viii) Centre of Approximated Interval (COAI) 

Centre of Area (COA) 

This is the most commonly used technique. This 
defuzzification can be expressed as 



v COA 



^^(x)^ 



where xcoaIs the crisp output, //^ (x) is the aggregated 

membership function and x is the output variable. 
This method is also known as center of gravity or 
centroid defuzzification method. The only disadvantage 
of this method is that it is computationally difficult for 
complex membership functions. 

Bisector of Area 

The bisector of area is the vertical line that divides the 
region into two sub-regions of equal area. The formula 
for xboa is given by 

x BOA a 4 

j M%( x ) dx = j Mg( x ) d x- 

a \ x BOA 

It is sometimes, but not always coincident with the 
centroid line. 

Largest of Maxima 

The maximum xwm takes the largest amongst all x that 
belong to [a 2 > a 3] • 

Smallest of Maxima 

The smallest output with the maximum membership 
function is taken as the crisp value and it is denoted by 

XSOM. 

Mean of Maxima 

In this method only active rules with the highest 
degree of fulfillment are taken into account. The 

output is computed as x M0M = ^(x L0M + x S0M ). 
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Regular Weighted Point (RWP) 

For the fuzzy number ft = (a^a^a^ , a 4 ) , the a -cut is 
A a = [L A (a),R A {a)~\ and the regular weighted point 
for A is given by (Saneifard 2012) 



rfM«)+M«r i 



J 



RWP(A) ■■ 



f [a)da 



]f(a]dt 



where /(«) = 



=j(L A (a) + R A (a)) f (a)da 
o 

\l-2a when a e [0,1/2] 



Via -1 when a e [1/2. , l] 

Graded Mean Integration Value of Fuzzy Number 

For the generalized fuzzy number A with 
membership function ju A (x) , according to Chen et al. 

(1999), the Graded Mean Integral Value P dGw (A) of A 
is given by 



M)- 



j" - w)L _1 (x) + wR' 1 (x)} dx 

o 

l 

j xdx 



l 

= 2j x{ (1 - hOZT 1 (x) + wR' 1 (x)| <& 



where the pre-assigned parameter w e [0, 1] refers to 
the degree of optimism, w = 1 represents an optimistic 
point of view, w = represents a pessimistic point of 
view and w = 0.5 indicates a moderately optimistic 
decision makers' point of view. 

Centre of the Approximated Interval (COAI) 

Let A be a fuzzy number with interval of confidence at 
the level a , then the a -cut is [A L (a),A R (a)] . The 

nearest interval approximation of A is with respect to the 



distance metric d is C d (A) : 
where 



| A L (a)da, j" A R (a)da 



d(A,B)= \\{A L {a)-B L {a)} 2 da + \{A R (a)-B R (a)} 2 da. 
V o o 
The interval approximation for the triangular fuzzy 

, ~ . \(a x +a 2 ) (a 2 +a 3 )l 
number A = (a l ,a 2 ,a i ) is — — , — ^ , for 

the parabolic fuzzy number A = (a 1 ,fl 2 ' a 3) is 

and for the trapezoidal fuzzy 



— (2^ +a 2 ), — (a 2 +2a 3 ) 



number A = (a x ,a 2 , a 3 ,a 4 ) is 



^■(fl! +a 2 ),^-(a 3 +a 4 ) 



The defuzzification formulas for different fuzzy 
numbers are given in the Table 1. Also, the positions of 
defuzzifying points of different defuzzification 
methods are shown in Fig. 2. 

TABLE 1 DEFUZZIFICATION FORMULAS FOR DIFFERENT FUZZY NUMBERS 



1 >.-lll//lli." 


1 IN 








I II N 






'j) 






lT . oa-oa 1 












BOA 




' 1 








MOM 












SOM 












LOM 


"a 










RWP 




«.) 








CMIV 




' 1 




+ 40,) 




luilhii 




COAI 




»,) 










^^OkT^ &MtV^^ 



FIG. 2 POSITION OF DEFUZZIFYING POINTS IN DIFFERENT 
DEFUZZIFICATION METHODS 

Mathematical Formulation of the Problem 

Let us consider a redundancy allocation problem of 
complicated reliability system with ten-component. 
Each subsystem is connected parallelly with identical 
components. Our objective is to maximize the overall 
system reliability subject to the given resource 
constraints. This can be done by finding the number of 
redundant components in each subsystem. 

The general form of the redundancy allocation 
problem in crisp form is as follows: 

Maximize R s (x) (1) 
subjected to g, (x) <b i ,i = l,2,--;tn , where 
x = (x l ,x 2 ,...,x n ) , 1 < L < Xj < iij, Xj is integer, /' = \,...,n , 
and b t is the z'-th available resource, i = 1, 2, • ■ •, m. 

Now, if the component reliabilities are imprecise, then 
the reliability of each subsystem and finally the overall 
system reliability will be imprecise. In this situation, 
the general form of the redundancy allocation problem 
can be written as follows: 

Maximize R s (x) (2) 
subjected to (x) < b { , i = 1, 2, ■ ■ ■, m where 
x = (x l ,x 2 ,...,x n ) , 1 < Ij < Xj < Uj, Xj is integer, j = \,...,n 

and b t is the i-th available resource which is imprecise , 
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i = l, 2,---,m. 

Fuzzy number representation of reliability of each 
component as well as different parameters of resource 
constraints is used to characterize the impreciseness of 
the system. 

Genetic Algorithm Based Constraints Handling 
Approach 

Evidently the problem (2) is a constrained optimization 
problem. Several techniques have been projected to 
deal with the constraints for solving the optimization 
problem. Recently, Gupta et al. (2009) and Bhunia et al. 
(2010) solved the optimization problem using Big-M 
penalty technique. In this method, the given 
constrained optimization problem is converted to an 
unconstrained one by penalizing a large positive 
number say, M . This method is known as Big-M 
penalty technique. In this work, the Big-M penalty 
technique has been used to handle the constraints. 

The corresponding unconstrained optimization problem 
is as follows: 

/ — \ defuz ( ) if x e S 
Maximize defuz I R s ) = \ y s ' (3) 
[-M if x <2 S 

where 

S = ^x:defuz(g j (x l ,x 2 ,...,x n )^ < defuz(b^,j = 1,2,... ,m j 

be the feasible space. Here, defuz^A^ denotes the 

defuzzified value of the fuzzy number A . The 
problem (3) is a non-linear unconstrained optimization 
problem with defuzzified objective of n integer 
variables In the next section, a solution 

procedure is prepared to solve the problem (3). 

Solution Procedure 

Since the problem in (3) is non-linear in nature with 
integer variables, it is difficult to reach an analytic 
solution (if any) to the problem (Gen and Cheng 1997). 
Furthermore, efficient treatment of integer non-linear 
optimization problem is one of the most difficult 
problems in practical optimization (El-Sharkawi 2008). 
Hence, to solve it, a meta-heuristic algorithm should 
be employed. 

A number of researchers have successfully used meta- 
heuristic algorithms to solve the complicated 
optimization problems in different fields of scientific 
and engineering disciplines. Some of these meta- 
heuristic algorithms are tabu search, genetic algorithm, 
particle swarm optimization, ant colony optimization 
and bee colony optimization among which genetic 



algorithm is the most popular one. 
Genetic Algorithm (GA) 

Genetic Algorithm (GA) is a stochastic search and 
optimization technique (Goldberg 1989). Gen and 
Cheng (1996) described the application of GA to 
combinatorial problems including reliability optimization 
problems. Perhaps, GA is the most widely known 
evolutionary computation method due to its simplicity, 
powerfulness and wide application. It works by the 
evolutionary principles and chromosomal processing 
in natural genetics. GA maintains a population P(f) 

for generation t, of chromosomes which are the set of 
genes, the part of solution. 

The most fundamental idea of Genetic Algorithm is to 
replicate the natural evolution process artificially in 
which populations undergo continuous changes 
through genetic operators, like crossover, mutation 
and selection. In particular, it is very useful to solve 
complicated optimization problems which cannot be 
solved easily by direct or gradient based mathematical 
techniques. It is very effective to handle large-scale, 
real-life, discrete and continuous optimization 
problems without making unrealistic assumptions and 
approximations. 

The approach of GA is applied to many engineering 
problem as well as decision making problems in 
various fields. Some distinguished characteristics of 
GA are as follows (Goldberg 1989): 

(i) GA works with a coding of solution set, not the 
solution themselves, 

(ii) GA searches over a population of solutions, not 
a single solution, 

(iii) GA uses payoff information, not derivatives or 
other auxiliary knowledge, 

(iv) GA applies stochastic transformation rules, but 
not deterministic. 

Each chromosome represents a potential solution to 
the problem under consideration and is evaluated to 
give some measure of its fitness. Some chromosomes 
undergo stochastic transformation by means of genetic 
operators to form new chromosomes. The genetic 
operator 'crossover' creates new chromosomes by 
combining some parts from two chromosomes and the 
other genetic operator 'mutation' produces offspring 
by making changes in a single chromosome. A new 
population is formed by selecting the chromosomes 
having better fitness value from the parent and the 
offspring populations taken together. After several 
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generations, the algorithm converges to the best 
chromosome having the best fitness value and it 
represents the optimal solution. 

The algorithm to implement GA is as follows: 

Step-1: Set population size (popsize), maximum 
number of generations (maxgen), probability of 
crossover (pcros), probability of mutation (pmute) and 
the bounds of decision variables. 
Step-2: Set t=0. 

Step-3: Initialize the chromosomes of the population 
P(t). 

Step-4: Compute the fitness function (objective 
function) for each chromosome of P (f) . 

Step-5: Find the chromosome having the best fitness 
value. 

Step-6: Set f=t +1. 

Step-7: If the termination condition is satisfied, then go 
to Step-13; otherwise the next step. 

Step-8: Select the population from the population 

P(t-1) of (t-1)- fh generation using tournament selection. 

Step-9: Apply crossover & mutation operators on P(f) 
to produce new population . 

Step-10: Compute the fitness function value of each 
chromosome of P[t ) . 

Step-1 1: Find the best chromosome from P(f) . 

Step-12: Find the better of the best chromosomes of 
P(t) & P(t-l) and store it; go to Step-6. 

Step-13: Print the best chromosome and its fitness 
value 

Step-14: End. 

The following vital components have been taken into 
account in the application of genetic algorithm. 

(i) GA parameters (population size, maximum 
number of generation, crossover rate and 
mutation rate) 

(ii) Chromosome representation 

(iii) Initialization of population 

(iv) Evaluation of fitness function 

(v) Selection process 

(vi) Genetic operators (crossover and mutation) 
GA Parameters 

In terms of implementing GA, the following GA 
parameters has been taken into account: 



Population size (popsize): Population size determines 
the amount of information stored by the GA. There is 
no clear rule how large it should be. The population 
size is problem dependent and will need to increase/ 
decrease with the dimension of the problem. 

Maximum number of generations (maxgen): It varies 
from problem to problem and depends upon the 
number of genes (variables) of a chromosome and it 
is prescribed to be the termination criterion for 
convergence of the solution. 

Probability of crossover (pcross): In Genetic Algorithms 
crossover is considered to be the main search operator. 
The crossover operator is used to thoroughly explore 
the search process. In crossover operator the genetic 
information between/among two or more individuals 
are blended to produce new individuals. Normal 
range of the crossover rate lies in [0.60, 0.95]. 

Probability of mutation (pmute): Mutation operator 
plays an important role in genetic algorithm. After 
crossover operation, mutation is performed. It is 
intended to prevent to the falling of all solutions in the 
population into a local optimum of the problem under 
consideration. Mutation operator randomly changes 
the offspring resulted from crossover. The mutation 
rate lies in [0.05, 0.20]. Sometimes mutation rate is 
considered to be 1/n where n is the number of genes 
(variables) of the chromosome. 

Chromosome Representation 

In the existing literature there are different types of 
representations of chromosomes, like, binary, real, 
octal, hexadecimal coding, amongst which real coding 
representations are well accepted. In this representation, 
a chromosome is coded in the form of vector/matrix of 
integer/floating point or combination of the both 
numbers and every component of that chromosome 
represents a decision variable of the problem. In this 
representation, each chromosome is encoded as a 
vector of integer numbers, with the same component 
as the vector of decision variables of the problem 
under consideration. This type of representation is 
accurate and more efficient as it is closed to the real 
design space and moreover, the string length of each 
chromosome is same as the number of decision 
variables. In this representation, for a given problem 
with n decision variables, a n-component vector 
x = (x l ,x 2 ,---,x n ) is used as a chromosome to represent 
a solution to the problem. A chromosome denoted as 
v k ( k = l,2,...,p_size) is an ordered list of n genes 
asv t = {v tl ,v t2 ,...,v b -,...,v fa }. 
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Initialization of Population 

After representation of chromosome, the next step is to 
initialize the chromosome that will take part in the 
artificial genetics. To initialize the population, first of 
all we have to find the independent variables and their 
bounds for the given problem. Then the initialization 
process produces population size number of 
chromosomes in which every component for each 
chromosome is randomly generated within the 
bounds of the corresponding decision variable. There 
are several procedures to select a random number of 
integer type. In this work, the following algorithm has 
been applied to select of an integer random number. 

An integer random number between a and b can be 
generated as either x = a + g or, x = b- g , where g is a 
random integer between 1 and \a -b\ . 

Evaluation of Fitness Function 

Evaluation of fitness function is same for natural 
evolution process in the biological and physical 
environments. After initialization of chromosomes of 
potential solutions, we need to see how relatively 
good they are. Therefore, the fitness value for each 
chromosome has to be calculated. In our work, the 
value of objective function of the reduced 
unconstrained optimization problems corresponding 
to the chromosome is considered as the fitness value of 
that chromosome. 

Selection of Fitness Function 

The selection operator which is the first operator in 
artificial genetics plays an interesting role in GA. This 
selection process is based on the Darwin's principle on 
natural evolution "survival of the fittest". The primary 
objective of this process is to select the above average 
individuals/chromosomes from the population 
according to the fitness value of each chromosome and 
eliminate the rest of the individuals/chromosomes. 
There are several methods to implement the selection 
process. 

In this work, the well known tournament selection 
with size two has been utilized. 

Genetic Operators 

After the selection process, other genetic operators, 
like crossover and mutation are applied to the 
resulting chromosomes those which have survived. 
Crossover is an operator that creates new individuals/ 
chromosomes (offspring) by combining the features of 
both parent solutions. It operates on two or more 



parent solutions at a time and produces offspring for 
next generation. In this work, intermediate crossover 
for integer variables has been used. 

The aim of mutation operator is to introduce the 
random variations into the population used to prevent 
the search process from converging to the local optima. 
This operator helps to regain the information lost in 
earlier generations and is responsible for fine tuning 
capabilities of the system and is applied to a single 
individual only. Usually, its rate is very low; because 
otherwise it would defeat the order building 
generated through the selection and crossover 
operations. In this work one-neighborhood mutation 
has been in use. 

Numerical Example 

The mathematical formulation of the redundancy 
allocation problem of a complex reliability system 
with ten subsystems is as follows 




FIG. 3 DIAGRAM OF A COMPLEX SYSTEM 
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TABLE 2 INPUT PARAMETERS VALUES FOR TFN AND PFN CASE 
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(4.3. 5.0. 
(4.3, 5.1. 
(5.5,6.8, 
(4.3, 5.4, 
(7.6. S.O. 



0.2 J 
7.h) 
0.3) 
8.7) 
7.8) 
6.2) 
(-. I i 
".4 I 
6.7) 

9.3) 



P - (207, 300,305), C = (306, 4(H), 4IO), W=(30K, 400. 404) 
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TABLE 3 INPUT PARAMETERS VALUES FOR TrFN CASE 



TABLE 4 RESULTS FOR TFN CASE 
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r 


n 
1 ■ 


1 


w 

■ 




(0.80.0.90,0.95. 0.98) 


14.1 1. 4.4,5.0. 5.5) 


16.9. 7.8. 8.5, 9.0) 


(7.5, 8.6, 9.2, 9.6) 


2 


(0.60, 0,70, 0,80, 0.8S) 


(4.5, 5.0,5.6,6.0) 


(5.4,6.3, 7.1,7.5) 


(5.3, 6.1,7.6, 8.3) 


3 


(0.85.0.90.0.%, 0.86) 


(3.0.3.2.4.0.4.5) 


(7.5, 8.1.9.0.9.21 


(7.2. 8.4. 9.3. 9.7) 


1 


(0.65.0.70, 0.SIJ, 0.86) 


(2.6,3.1,18.4.3) 


(8.0.9.2. 9.7,9.8) 


(6.4. 7.5,8.7, 9.1) 


5 


(0.62.0.70,0.79,0.82) 


(3.8, 4.4,5.3.6.0) 


(6.6. 7.3. 8.3, 9.01 


(5.4. 6.7, 7.8, 8.6) 


6 


(0.75.0.90,0.95,0.97) 


(2.1,2.9,3.4.4.0) 


(5.7. 6.2. 7.7, 83) 


(4.3, 5.0, 6,2, 7.4) 


7 


(0.50.0.70.0.79, 0.83) 


(2.4. 2.9.3.7.4.2) 


15.6.6.5. 7.3.8.0) 


(4.3. 5.1.6.4.7.5) 


8 


(0.76.0.80.0.87, 0.91) 


(2.6.3.1.4.3.4.6) 


(7.3. 8.1.8.7.8.9) 


(5.5. 6.8.7.4.8.3) 


9 


(0.75. 0.80, 0.S6, (190) 


(1.8, 2.5, 3.5. 3.8) 


(S.5,9.2. 10.0, 102) 


(4.3. 54, 6.7, 7.2) 


in 


(0.70.0.80,0.95,0.97) 


(1.7, 2.8.34.40) 


(6.0. 7.9. 8.5, 9.0) 


(7.6. 8.9, 9.. 1 9.8) 


P = 1 297. 300,305308), C = (3%. 400. 410,413), H'=(39S. 400. 404,407) 



In all the cases, the problem has been fixed by real 
coded advanced genetic algorithm. In this algorithm, 
we have used tournament selection, intermediate 
crossover and one neighborhood mutation as genetic 
operators. For this purpose, the code has been 
prepared for this algorithm in C Programming 
language. The corresponding computational work has 
been done on a PC with Intel Core-2 duo processor in 
LINUX environment. For each problem, 20 
independent runs have been performed to determine 
the best found system reliability which is nothing but 
the optimal value of system reliability. In this 
computation, the values of genetic parameters, like 
popsize, maxgen, pcross and pmute have been taken as 
300, 300, 0.85 and 0.15 respectively. 

The computational results have been shown in Tables 
4-6 for different parametric values. From Table 4, it 
follows that in the case of TFN parametric values, the 
best optimum system reliability is obtained in the 
cases of MOM, SOM and LOM. In case of PFN 
parametric values, the best optimum system reliability 
is obtained in cases of MOM, SOM and LOM as in the 
case of TFN parametric values and given in Table 5. 
Table 6 represents the results in the case of TrFN 
parametric values and indicates that the best system 
reliability is available in the case LOM. All the results 
have been shown in Fig. 4. 
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TABLE 5 RESULTS FOR PFN CASE 
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FIG. 4 SYSTEM RELIABILITY FOR DIFFERENT 
DEFUZZIFICATION METHODS 



Concluding Remarks 

The reliability of a component in a system might be 
imprecise rather than precise for several reasons. This 
impreciseness may be represented by diverse ways. In 
this paper, we have represented this by fuzzy number. 
Then the problem has been converted into crisp 
nonlinear programming problem after defuzzification 
in which the objective functions and all the constraints 
are crisp valued. For further research, one may apply 
other heuristic methods, like, DE, SA, PSO, etc. to 
solve the problem discussed in this paper. Also, the 
proposed technique may be applied in solving the 
realistic engineering and other optimization problems. 
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